APPENDICIES 



Appendix 1: Abbreviation Key for Mass- Volume Model 



Abbreviation 



Kf sd = associated rate constant for stomach and duodenum 



Ka dj = associated rate constant for duodenum and jejunum 



Ka ji = associated rate constant for jejunum and ileum 



Ka ie = associated rate constant for ileum and colon 



Ka co = associated rate constant for colon and excretion 



SD trans = transfer rate between stomach and duodenum 



DJ trans = transfer rate between duodenum and jejunum 



JL trans = transfer rate between jejunum and ileum 



IC trans = transfer rate between ileum and colon 



Waste = transfer rate between colon and excretion 



pH s = pH stomach 



pH s2 = pH duodenum 



pH s3 = pH jejunum 



pH s4 = pH ileum 



pH s5 = pH colon 



sol profile = solubility profile for stomach 
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sol profile 2 = solubility profile for duodenum 



sol profile 3 = solubility profile for jejunum 



sol profile 4 = solubility profile for ileum 



sol profile 5 = solubility profile for colon 



stom ka = associated rate constant for stomach compartments 1 and 2 



duo ka = associated rate constant for duodenum compartments 1 and 2 



Jej ka = associated rate constant for jejunum compartments 1 and 2 



II ka = associated rate constant for ileum compartments 1 and 2 



Colon ka = associated rate constant for colon compartments 1 and 2 



S A stom = surface area of stomach 



SA duo = surface area of duodenum 



SA jej = surface area of jejunum 



S A il = surface area of ileum 



S A colon = surface area of colon 



Perm stom = permeability of stomach 



Perm duo = permeability of duodenum 



Perm jej = permeability of jejunum 



Perm il = permeability of ileum 



Perm colon = permeability of colon 
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Ka sd = associated rate construct for stomach fluid absorption 

Ka du = associated rate construct for duodeunm fluid absorption 

Ka je = associated rate construct for jejunm fluid absorption 

Ka il = associated rate construct for ileunm fluid absorption 

Ka co = associated rate construct for colon fluid absorption 

Note: other abbreviations adhere to above descriptors and are self 
explanatory 
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Appendix 2: Equations, Parameters and Values For Mass-Volume Model 

amt_plasma(t) = amt_plasma(t - dt) + (trans_21 + ka - elimination - trans_12) * dt 
IN IT amt_plasma = 0 

INFLOWS: 

trans_2 1 = k21*comp_2 
ka = tot_abs_rate 

OUTFLOWS: 

elimination = amt_plasma* k_elim 
trans_12 = kl2*amt_plasma 

blood_side_col(t) = blood_side_col(t - dt) + (colon_ka_5) * dt 
INIT blood_side_col = 0 

INFLOWS: 

colon_ka_5 = IF Vol_colon*sol_profile_5 >=Colon THEN Colon*SA_colon*perm_colon*3600 
ELSE Vol_colon*sol_profile_5*SA_colon*perm_colon*3600 
blood_side_dou(t) = blood_side_dou(t - dt) + (duo_ka) * dt 
INIT blood_side_dou = 0 

INFLOWS: 

duo_ka = IF Vol_duod*sol_profile_2 >= duodenum THEN 
duodenum*SA_duo*perm_duo*3600 ELSE Vol_duod* sol_profile_2* S A_duo*perm_duo* 3 600 
blood_side_il(t) = blood_side_il(t - dt) + (Il_ka) * dt 
INIT blood_side_il = 0 

INFLOWS: 

Il_ka = IF Vol_ileum*sol_profile_4 >=Ueum THEN Ileum*SA_Il*perm_Il*3600 ELSE 
Vol_ileum*sol_profile_4*SA_Il*perm_Il*3600 
blood_sideJej(t) = blood_sideJej(t - dt) + (Jej_ka) * dt 
INIT blood_sideJej = 0 

INFLOWS: 

Jej_ka = IF VolJej*sol_profile_3 >=Jejunum THEN Jejunum*SAJej*permJej *3600 ELSE 
Vol Jej *sol_profile_3 * SA Jej *perm Jej *3600 
blood_side_sto(t) = blood_side_sto(t - dt) + (stom_ka) * dt 
INIT blood_side_sto - 0 

INFLOWS: 

stom_ka = IF Vol_stom*sol_profile >= Stomach THEN Stomach*SA_stom*perm_stom*3600 
ELSE Vol_stom*sol_profile*SA_stom*perm_stom*3600 
Colon(t) = Colon(t - dt) + (IC_trans - Waste - colon_ka_5) * dt 
INIT Colon = 0 
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INFLOWS: 

ICtrans = ka_ic*Ileum 

OITFLOWS: 

Waste = ka_col* Colon 

coion_ka_5 = IF Vol_colon*sol_profile_5 >=Colon THEN Colon*SA_colon*perm_colon*3600 
ELSE Vol_colon*sol_profile_5 * SA_colon*perm_colon*3600 
comp_2(t) = comp_2(t - dt) + (trans_12 - trans_21) * dt 
1NIT comp_2 = 0 

INFLOWS: 

trans_12 = kl2*amt_plasma 

OUTFLOWS: 

trans_21 = k21*comp_2 

duodenum(t) = duodenum(t - dt) + (SDjrans - duo_ka - DJ_trans) * dt 
INIT duodenum = 0 

INFLOWS: 

SDjrans = if Stomach >0 then kf_sd* Stomach else 0 

OUTFLOWS: 

duo_ka = IF VoI_duod*sol_profile_2 >= duodenum THEN 

duodenum*SA_duo*perm_duo*3600 ELSE Vol_duod*sol_profile_2*SA_duo*perm_duo*3600 

DJtrans = ka_dj*duodenum 

excretion(t) = excretion(t - dt) + (vol_cw) * dt 

INIT excretion = 0 

INFLOWS: 

vol_cw = Vol_colon*ka_col 

excretion_2(t) = excretion_2(t - dt) + (Waste) * dt 

INIT excretion_2 - 0 

INFLOWS: 

Waste = ka_col*Colon 

Ileum(t) = Ileum(t - dt) + (JLjrans - ICjrans - Il_ka) * dt 
INIT Ileum = 0 

INFLOWS: 

JL_trans = kaji* Jejunum 

OUTFLOWS: 
ICtrans = ka_ic* Ileum 

Il_ka = IF Vol_ileum*sol_profile_4 >=Ileum THEN Ileum*SA_Il*perm_Il*3600 ELSE 

Vol_iIeum*sol_profile_4*SA_Il*perm_Il*3600 

Jejunum(t) = Jejunum(t - dt) + (DJjrans - JL_trans - Jej_ka) * dt 
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INIT Jejunum = 0 



INFLOWS: 

DJ_trans = ka_dj*duodenum 

OUTFLOWS: 

JLjrans = kaJi*Jejunum 

Jej_ka = IF VolJej*sol_profile_3 >=Jejunum THEN Jejunum*SAJej*permJej *3600 ELSE 

VolJej*sol_profile_3*SAJej*permJej*3600 

serosal_col(t) = serosal_col(t - dt) + (Adsorp_col - col_secretion) * dt 

INIT serosal_col - 0 

INFLOWS: 

Adsorp_col = PULSE(1.67,0,.l)+0*Vol_colon*ka_co 

OUTFLOWS: 
colsecretion = 0 

serosal_dou(t) = serosal_dou(t - dt) + (Adsorp_Duo - duo_secretion) * dt 
INIT serosal_dou = 0 

INFLOWS: 

Adsorp_Duo = PULSE(1 0.82,0,. l)+0*Vol_duod*ka_du 
OUTFLOWS: 

duo_secretion = PULSE(1 0.82,0,. 1) 

serosal_ill(t) = serosal_ill(t - dt) + (Adsorpt_ill - ile_secretion) * dt 
INIT serosal_ill = 0 

INFLOWS: 

AdsorpUll = PULSE(8.83,0,.10)+0*Vol_ileum*ka_il 
OUTFLOWS: 

ile_secretion = PULSE(1.50,0,.l) 

serosal Jej(t) = serosal Jej(t - dt) + (Adsorpjej - jej_secretion) * dt 
INIT serosal J ej = 0 

INFLOWS: 

Adsorpjej = PULSE(1 5.76,0,. l)+0*VolJej*kaJe 
OUTFLOWS: 

jej_secretion = PULSE(2.67,0,.l) 

serosal_sto(t) = serosal_sto(t - dt) + (Adsorp_Stom - Stom_Secretion) * dt 
INIT serosalsto = 0 

INFLOWS: 

AdsorpStom = 0*Vol_stom*ka_sd 
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OUTFLOWS: 

Stom_Secretion = PULSE(1 6.67,0,- 1) 

Stomach(t) = Stomach(t - dt) + (- SDjrans - stom_ka) * dt 

INIT Stomach =1000 

OUTFLOWS: 

SDtrans = if Stomach >0 then kf_sd* Stomach else 0 

stom ka = IF Vol_stom*sol_profile >= Stomach THEN Stomach*SA_stom*perm_stom*3600 
ELSE Vol_stom*sol_profile*SA_stom*perai_stom*3600 
total_drug_absorbed(t) = total_drug_absorbed(t - dt) + (tot_abs_rate) * dt 
INIT totaI_drug_absorbed = 0 

INFLOWS: 

tot_abs_rate = stom_ka+duo_ka+Jej_ka+Il_ka+colon_ka_5 
Total_Elimination(t) = Total_Elimination(t - dt) + (elimination) * dt 
INIT Total_Elimination = 0 

INFLOWS: 

elimination = amt_plasma*k_elim 

Vol_colon(t) = Vol_colon(t - dt) + (vol_ij + col_secretion - vol_cw - Adsorp_col) * dt 
INIT Vol_colon = 0 

INFLOWS: 

voljj = Vol_ileum*ka_ic 
col_secretion = 0 

OUTFLOWS: 

vol_cw = Vol_colon*ka_col 

Adsorp_col = PULSE(1.67,0,.l)+0*Vol_colon*ka_co 

Vol_duod(t) = Vol_duod(t - dt) + (vol_sd + duo_secretion - voil_dj - Adsorp_Duo) * dt 
INITVol_duod = 0 

INFLOWS: 

volsd = kf_sd*Vol_stom 
duo_secretion = PULSE(10.82,0,.l) 

OUTFLOWS: 

voil_dj = Vol_duod*ka_dj 

Adsoip_Duo = PULSE(10.82,0,.l)+0*Vol_duod*ka_du 

Vol_ileum(t) = Vol_ileum(t - dt) + (volji + ile_secretion - Adsorpt_ill - vol_ij) * dt 
INIT Vol_ileum = 0 

INFLOWS: 

volji = VolJej*kaJi 

ile_secretion = PULSE(1 .50,0,.!) 
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OUTFLOWS: 

Adsorpt_ill = PULSE(8.83,0 v 10)+0*Vol_ileum*kaJl 
voljj = Vol_ileum*ka_ic 

Voljej(t) = VolJej(t - dt) + (voil_dj + jej_secretion - volji - Adsorpjej) * dt 
INIT Voljej = 0 

INFLOWS: 

voil_dj = Vol_duod*ka_dj 
jej_secretion = PULSE(2.67,0,.l) 

OUTFLOWS: 
volji = VolJej*kaJi 

Adsorpjej = PULSE(15.76 J 0 v l)+0*VolJej*kaJe 

Vol_stom(t) = Vol_stom(t - dt) + (Stom_Secretion - vol_sd - Adsorp_Stom) * dt 
INIT VoLstom = PULSE(8.33,0 V 1) 

INFLOWS: 

Stom_Secretion = PULSE(16.67,0 V 1) 

OUTFLOWS: 

vol_sd = kfsd*Vol_stom 

Adsorp_Stom = 0* Vol_stom*ka_sd 

concjplasma = (amt_plasma/volurne)* mg_to_ug 

kl2 = .839 

k21 = .67 

ka_co = 1 

ka_col = 3 

ka_dj = 3 

ka_du = 1 

ka_Jc = 3 

ka_il = 8.83 

ka je = 1 

kaji = 3 

ka_sd = 1 

kf_sd = 2.8 

k_elim = .161 

mg_to_ug= 1000 

perm_colon = 3.80e-6 

perm_duo= 1.10e-6 

perm_Il = 4.06e-6 

permjej =2.17e-6 

perm_stom = 1 .10e-6 

ph_s=1.5 

ph_s_2 = 6.6 

ph_s_3 = 6.6 
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ph_s_4 = 7.5 
ph_s_5 = 6.6 
SA_colon= 138 
SA_duo = 125 
SAJ1 = 102 
SAJej = 182 
SA_stom = 50 
volume = 4*19200 
sol_profile = GRAPH(ph_s) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4.50, 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65), 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_2 = GRAPH(ph_s_2) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00. 4.00), (3.50, 3.80), (4.00, 3.65), (4.50, 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65), 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_3 = GRAPH(ph_s_3) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4.50, 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65), 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_4 * GRAPH(ph_s_4) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4.50. 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65), 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
sol_profile_5 = GRAPH(ph_s_5) 

(1.00, 63.0), (1.50, 25.0), (2.00, 10.0), (2.50, 5.00), (3.00, 4.00), (3.50, 3.80), (4.00, 3.65), (4.50, 
3.50), (5.00, 3.65), (5.50, 3.65), (6.00, 3.65), (6.50, 3.65), (7.00, 3.65), (7.50, 3.65), (8.00, 3.65), 
(8.50, 4.00), (9.00, 5.00), (9.50, 12.0), (10.0, 23.5) 
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Appendix 3: Abbreviation Key For GI Model 



The legend/key has been divided into sub-sections corresponding to the sub-sections of the 
model diagram. 

Numbered suffixes (1,2, 3, 4, 5, 6) have been assigned to distinguish between intestinal regions 
(stomach, duodenum, jejunum, ileum, colon, and waste, respectively). 

1 - stomach 

2 - duodenum 

3 - jejunum 

4 - ileum 

5 - colon 

6 - waste 

For example, VOL 1 is the volume in the stomach, MASS 3 is the insoluble mass in the jejunum. 
In the equations, COMP 1 indicates the stomach, COMP 2 the duodenum, COMP 3, the jejunum, 
etc. 

Ghosts are listed under the sub-section containing the original reservoir, flow regulator, or 
converter. 

Abbreviations listed in italics are regionally dependent and set up as arrays to allow independent 
values for each intestinal region. 

In general, ADJ as a prefix indicates a calculated parameter value (ADJ = adjusted), while ADJ 
as a suffix indicates an adjustment parameter (ADJ = adjustment). 



Intestinal model 

Reservoirs/Compartments 

Fluid volume absorbed 
Fluid volume 

Mass of drug contained with a formulation or controlled release 
device 

Insoluble mass of drug (not contained within the formulation or 
controlled release device) 
Soluble mass of drug 
Mass of drug absorbed 

Flow regulators 



VOL ABS 

VOL 

CREL 

MASS 

SOL 

ABSORPTION 
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RE AB S Rate of water absorption 

VOL OUT Fluid volume transit rate 

CR OUT Formulation or controlled release device transit rate 

CR INPUT Drug release rate from formulation or controlled release device 

MASS OUT Insoluble drug mass transit rate 

DISS PRECIP Dissolution rate 

SOL OUT Soluble drug mass transit rate 

FLUX Absorption rate 



AD J PARMS (Adjustment Parameters) 

VOL AD J Fluid volume absorption adjustment parameter 

DISS ADJ Dissolution rate adjustment parameter 

TRANSIT ADJ Transit time adjustment parameter 

S A ADJ Surface area adjustment parameter 

FLUX ADJ Passive Absorption adjustment parameter 

EFFLUX ADJ Efflux or secretion adjustment parameter 

CARRIER ADJ Active absorption adjustment parameter 



PARMS (Parameters) 



VOL PARM 
SURFACE AREA 
DOSE 

INIT VOLUME 
TIME IN HOURS 
pH 

PARACELLULAR 



Fluid volume absorption rate constant 
Surface area available for absorption 
The administered dose of drug 
The administered volume of water or fluid 
A clock 

The physiological pH value 

A user controlled switch used to adjust absorption based on 
absorption mechanism 



TRANSIT TIME 



TRANSFERS 
CUMUTT 

ADJ TRANSIT TIME 
USER TT INPUT 



GI transit rate constant 
Cumulative transit time 

Adjusted GI transit time incorporating adjustment parameter and 
user input 

User controlled adjustments to the GI transit time 



OUTPUT CALCULATIONS 

ABSORBED TOTAL Total mass of drug absorbed (sum of ABSORPTION 1...5) 
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FD P 0/o Fraction or the dose absorbed into portal vein x 1 00 

FLUX TOTAL Total absorption rate (sum of FLUX 1 ...5) 

CUM DISS Cumulative drug mass dissolved 

CR Release Cumulative drug mass released from formulation 

CUM DISS RATE Sum of DISS PRECIP 1...5 

CR cumrate Summ of CR INPUT 1 ...5 



PERMEABILITY CALCULATION 



ADJ PERM 

ACT PE 
Km 

REGIONAL 

PASS PE 
RC 

RCSUM 



Adjusted permeability ncorporating all transport mechanisms and 

relevant adjustment parameters 

Active or carrier-mediated absorptive permeability 

Constant from the Michaelis-Mentin type permeability equation for 

active transport 

Passive permeability after regional correlation calculation (same as 
PASS PE if regional correlation is not used) 
Passive permeability entered by user 

A logical function used in determining the regional correlation ' 
A logical function used in determining the regional correlation 



SOLUBILITY CALCULATION 

User supplied pH value for which a solubility value is available 
User supplied solubility value corresponding to the USER pH value 
Solubility calculated (if necessary) at the appropriate pH value 
using the entered USER pH and USER SOLUB values 



CONTROLLED RELEASE CALCULATION 

The instantaneous release rate from the formulation 
The total dose contained with the formulation 
The cumulative drug mass release profile 
The cumulative drug mass release profile 

Note: CR AT TIME holds the value at the current time value (t), CR AT LAST holds the value 
at the immediately proceeding time value (t-1) 



USERpH 
USER SOLUB 
ADJ SOLUB 



CR RATE 
CRDOSE 
CR AT TIME 
CRAT LAST 



CONC CALCULATION 

CONCENTRATIONS The dissolved drug concentration 
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DISSOLUTION CALCULATION 

PRECIP Precipitation rate constant 

DISSOL Dissolution rate constant 

ADJ DISS PRECIP Adjusted rate constant incorporating PRECIP, DISSOL and 

calculated concentration 



O 
m 

m 
m 
m 

m 

w 

y 
ru 

Q 
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Appendix 4: Equations, Parameters and Values For GI Model 



AD J PARMS 

O CARRIER_ADJ[COMPS] = 0 
3 0ISS_ADJ(COMP_1] = 1 
Q DISS_ADJ[COMP_2] = 1 
Q DISS_ADJ[COMP_3] = 1 
3 DISS_ADJ[COMP_4] = 1 
Q DISS_ADJ[C0MP_5] = 1 
Q EFFLUX_ADJ[COMPS] = 1 
Q FLUX_ADJ[COMPJls1 
O FLUX_ADJ[COMP_2] = 10 
0 FLUX_ADJ[COMP_3] = 8 
0 FLUX_ADJ[COMP_4] = 2 
O FLUX_ADJ[COMP_5] = 1 
Q SA_ADJ[COMP_1] = 1 
0 SA_ADJ[COMP_2] = 1 
O SA_ADJ{COMP_3] = 1 
O SA_ADJ[C0MP_4] = 1 
0 SA_ADJlCOMP_5] * 1 
0 TRANSIT_ADJ[COMPJ] = 1 
O TRANSIT_ADJ[COMP_2] = 1 
Q TRANSIT_ADJ[COMP_3] « 1 
O TRANSIT_ADJ[COMP_4] * 1 
O TRANSIT_ADJ[COMP_5] « 1 
Q VOL_ADJ[COMP_1] = 1 
0 VOL_ADJ[COMP_2] = 1 
0 VOL_ADJ[COMP_3] = 1 
O VOL_ADJ[COMP_4] = 1 
O VOL_ADJ{COMP_5] = 1 
O CONC CALCULATION 

O CONCENTRATIONS[COMP_1] = if VOL_1 =0.0 then 0 else if 

ADJ SOLUB(COMP_1]<SOL_1/VOL_1 then ADJ_SOLUB[COMP_1] else SOL_1A/OL_1 + 

0*(SOL 2+SOL_5+SOL_3+SOL_4+VOL_3+VOL_2+VOL_4+VOL_5) 
O CONCENTRATIONS[COMP_2] * if VOL_2 = 0.0 then 0 else if (VOL_2<1e-6 AND SOL_2<1e-7) 

then 0 else if ADJ_SOLUB[COMP_2]<SOL_2A/OL_2 then ADJ_SOLUB[COMP_2] else 

SOL_2/VOL_2 

♦0*(SOL_1+SOL 5+SOL 3+SOL_4+VOL_3+VOL_1+VOL_5+VOL_4) 
O CONCENTRATIONS[COMP_3] =lf VOL_3 = 0.0 then 0 else if (VOL_3<1e-6 AND SOL_3<1e-7) 
then 0 else if ADJ_S0LUB[C0MP_3]<S0LJW0L_3 then ADJ_SOLUB[COMP_3] else 
SOL_3A/OL_3 

+0*(SOL_1 +SOL_2+SOL 4+SOL_5+VOL_5+VOL_4+VOL_1 +VOL_2) 
O CONCENTRATIONS[COMP_4] -If VOL_4 = 0.0 then 0 else if (VOL_4<1 e-6 AND S0L_4<1 e-7) 
then 0 else if ADJ_SOLUB[COMP_4]<SOL_4A/OL_4 then ADJ_SOLUBfCOMP_4] else 
SOL_4/VOL_4 

+0*(SOL J +SOL_2+SOL_3+SOL_5+VOL_1 +VOL_2+VOL_3+VOL_5) 



lib 



" C0NCENTRATIONS{COMP_5] * if V0L_5 = 0.0 then 0 else if (VOL_5<1e-6 AND SOL_5<le-7) 
then 0 else if ADJ_SOlUB[COMP_5}<SOL_5/VOL_5 then ADJ_SOLUB(COMP_5] else 
SOL_5/VOL_5 

♦0*(SOL_1 +SOL_4+SOL_3+SOL_2+VOL_3+VOL_1 +VOL_2+VOL_4) 
CONTROL RELEASE CALCULATION 
Q CR_DOSE = 0 

r CR_RATE = (CR_AT_TIME-CR_AT_LAST)*20*(CR_DOSE/100) 

?» CR_AT_LAST = GRAPH(TIME-DT) 

(0.00, 0.00), (0.25, 17.7), (0.5, 31.5), (0.75, 42.2), (1.00, 50.6). (1.25, 57.1), (1.50, 62.1), (1.75, 
66 1), (2.00, 69.2), (2.25, 71.6), (2.50. 73.4), (2.75, 74.9), (3.00, 76.0). (3.25, 76.9), (3.50, 77.6), 
(3.75. 78.1), (4.00. 78.5), (4.25. 78.9), (4.50. 79.1), (4.75. 79.3), (5.00, 79.5). (5.25, 79.6), (5.50, 
79.7), (5.75. 79.7), (6.00, 79.8), (6.25, 79.8). (6.50, 79.9). (6.75, 79.9), (7.00, 79.9), (7.25, 79.9), 
(7.50, 80.0), (7.75, 80.0), (8.00. 80.0), (8.25, 80.0), (8.50. 80.0). (8.75, 80.0). (9.00, 80.0), (9.25. 
80 0), (9.50, 80.0), (9.75, 80.0), (10.0, 80.0), (10.3, 80.0), (10.5. 80.0), (10.8, 80.0), (11.0. 80.0). 
(11.3. 80.0). (11.5, 80.0), (11.8, 80.0). (12.0. 80.0), (12.3, 80.0). (12.5, 80.0), (12.8, 80.0), (13.0, 

80.0) ... 

0 CR AT TIME = GRAPH(TIME) 

(0.6b, 0~00). (0.25. 17.7), (0.5, 31.5). (0.75. 42.2). (1.00. 50.6), (1.25, 57.1), (1.50, 62.1), (1.75. 

66.1) , (2.00, 69.2), (2.25, 71.6), (2.50. 73.4). (2.75. 74.9). (3.00. 76.0). (3.25. 76.9), (3.50. 77.6), 
(3.75, 78.1), (4.00, 78.5), (4.25. 78.9). (4.50, 79.1). (4.75. 79.3). (5.00. 79.5). (5.25. 79.6). (5.50, 
79.7), (5.75. 79.7). (6.00, 79.8), (6.25. 79.8). (6.50, 79.9), (6.75, 79.9), (7.00. 79.9). (7.25. 79.9). 
(7 50, 80.0). (7.75, 80.0), (8.00, 80.0). (8.25, 80.0), (8.50, 80.0). (8.75, 80.0). (9.00, 80.0), (9.25, 
80.0), (9.50, 80.0), (9.75, 80.0). (10.0, 80.0), (10.3,80.0), (10.5. 80.0). (10.8, 80.0), (11.0, 80.0), 
(1 1.3, 80.0), (11.5, 80.0). (11.8. 80.0), (12.0, 80.0). (12.3. 80.0). (12.5, 80.0), (12.8. 80.0), (13.0, 
80.0)... 

DISSOLUTION CALCULATION 

O ADJ_DlSS_PRECIP[COMP_1] * if VOL J *0 then 0 else if 
(SOL_1 /v6L_1 <AD J_SOLUB(COMP_lf) then 

(DISSOL[COMP_1]*DISS ADJ[COMPJ] # MASS J*(ADJ_SOLUB[COMP_1]-SOL_1/VOLJ )) else 
((SOL 1/V0L_1>ADJ_S6LUB[C0MP_1])*PRECTP[C0MPJ]+ 

0*(MASS_1+MASS_2+MASS_3+MASS_4+MASS_5+SOL_1+SOL_2+SOL_3+SOL_4+SOL_5W 
OL 1+VOL 2+VOL_3+VOL 4+VOL.5) 
O ADJ_DISS_PRECIPlCOMP_2] = if VOL_2=0 then 0 else if 
(SOL_2/VOL_2<ADJ S0LU~BIC0MP_21) then 

(DISSOL(COMP 2rDISS_ADJ[COMP_2] , MASS_2*(ADJ_SOLUB[COMP_2)-SOL - 2A/OL_2))else 
((SOL 2/VOL 2VADJ SOLUBfCOMP 2])*PREc7pfCOMP_21 

+0'(MASS_1+N^SSj+MASSJ+MASS_4+MASS.5+SOL_1+SOL_2+SOL_3+SOL_4+SOL_5+V 
0L_1+V0C_2+v6l_3 + V0L_4+V0L_5) ~ 
O ADJ DiSSjPRECIP[COMP_33 = ifVOL_3=0 then 0 else if 
(SOL_3/VOL 3<ADJ SOLUBlCOMP_3]j then 

(DISSOL[COMP 3]*DISS ADJ[COMP 3rMASS_3*(ADJ_SOLUB[COMP_3]-SOL_3A/OL_3)) else 
((SOL_3/VOL 3>ADJ SOLUB[COMP 3])*PRECIP{COMP_3] 

+0 , (MASS_1^MASS_2+MASS_3+MASS_4+MASS_5+SOL_1+SOL_2+SOL_3+SOL - 4*SOL_5+V 
OLJ+VOL 2+VOL 3+VOL_4+VOL_5) ~ 
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"~ ADJ_DISS_PREC1P{C0MP_4] = if VOL_4 S 0 then 0 else if 
"~" iSOL_4/VOL_4<ADJ_SOLUB[COMP_4]) then 

(DISSbL[COMP_4]*DISS_ADJ[COMP_4fMASS_4*(ADJ_SOLUB[COMP_4]-SOL_4A/OL 4)) else 

((S0L_4A/0L_4^ADJ - S0LUB[C0MP_4])*PRECTP{C0MP_4] 

+0*(MASSJ +MASS_2+MASS_3+MAis_4+MASS_5+SOL_1 +S0L_2+S0L 3+SOL_4+SOL 5+V 
0L_1 +VOL_2+VOL_3+VOL_4+VOL_5) 
C ADJ_DISS_PREClPtC0MP_5] = if VOL_5=0 then 0 else if 
(SOL_5A/OL_5<ADJ_SOLUB[COMP_5]) then 

(DlSSOL[COMP 5]*DlSS_ADJ[COMP_5rMASS_5*(ADJ_SOLUB[COMP_5]-SOL_5A/OL_5))else 
({SOL_5/VOL_5£ADJ_SOLUB[COMP_5]rPRECIP[COMP_5] 

+0'(MASS_1+MASS_2+MASS_3+MASS - 4+MASS_5+SOL,1+SOL_2+SOL_3+SOL_4+SOL_5+V 

OLJ+VOL 2+VOL_3+VOL_4+VOL_5) 
O DISSOLICOMPJ] = 1 
O DISSOL[COMP_2] = 1 
r DISSOL[COMP_3] - 1 
■5 DISSOL[COMP_4] = 1 
5 DISSOL{COMP_5] = 1 
C PRECIP[COMPJ] = 10 
Q PRECIP[COMP_2] = 10 
3 PRECIP[COMP_3] = 10 
C PRECIPfCOMPj*] = 10 
£ PRECIP[COMP_5] = 10 
INPUTS 

INTESTINAL MODEL 

□ ABSORPTIONJ (t) = ABSORPTIONJ (t - dt) + (FLUXJ ) • dt 
[NIT ABSORPTIONJ = 0 

INFLOWS: 
■*> FLUX 1 = 

CONCENTRATIONS[COMP_1]*ADJ_PERM[COMP_irSURFACE_AREA[COMP_1] 

□ ABSORPTION_2(t) = ABSORPTION_2(t - dt) + (FLUX.2) • dt 
INIT ABSORPTION_2 = 0 

INFLOWS: 

FLUX 2 = 

CONCENTRATlONS(COMP_2rADJ_PERM(COMP_2rSURFACE_AREA[COMP_2] 

□ ABSORPTION_3(t) = ABSORPTION_3(t - dt) + (FLUX.3) * dt 
INITABSORPTIQN_3*0 

INFLOWS: 

FLUX 3 = 

CONCENTRATIONS[COMP_3rADJ_PERM[COMP_3rSURFACE_AREA[COMP_3] 

□ ABSORPTION_4(t) = ABSORPTION_4(t - dt) + (FLUX.4) * dt 
INIT ABSORPTIONS * 0 

INFLOWS: 
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# FLUX_4 = 

CONCENTRAT]ONS[COMP_4]*ADJ_PERM[COMP_4]«SURFACE_AREA[COMP_4] 

d ABS0RPTI0N_5(t) = ABS0RPTI0N_5(t - dt) + (FLUX_5) * dt 
INIT ABS0RPTI0N_5 = 0 
INFLOWS: 
^ FLUX_5 = iftime<32then 

CONCENTRATIONS[COMP_5rADJ_PERM[COMP_5rSURFACE_AREA(COMP 5] l (32-ti 
me)/48'(VOL_5/1 7.2) else 0 
C REL 1(t) = C_REL_1(t - dt) ♦ (- CR OUT 1 - CRJNPUTJ ) * dt 
INIT C_REL_1 = CR_DOSE 
OUTFLOWS: 

CR_OUT_1 = IF TIME >= CUMU_TT[COMP_1] THEN C_RELJ*10000 ELSE 0 

# CR_INPUT_1 = if TIME>CUMUJT(COMP_1] then 0 else CR_RATE 
C C_REL_2(t) = C_REL_2(t - dt) ♦ (CR_OUT_1 - CR_OUT_2 - CRJNPUT.2) * dt 

INIT C_REL_2 = 0 
INFLOWS: 

CR_OUT_1 = IF TIME >= CUMUJTICOMPJ] THEN C_RELJ*10000 ELSE 0 
OUTFLOWS: 

CR_OUT_2 = IF TIME >= CUMU_TT[COMP_2] THEN C_REL_2*10000 ELSE 0 
CR_INPUT_2 = if TIME>CUMU_TT[COMP_2] then 0 else CR_RATE 

□ C_REL_3(t) = C_REL_3(t - dt) + (CR_OUT_2 • CrToUT_3 • CRJNPUT 3) * dt 
INIT C_REL_3 = 0 

INFLOWS: 

CR_OUT_2 = IF TIME >= CUMU_TT[COMP_2] THEN C_REL_2*10000 ELSE 0 
OUTFLOWS: 

^ CR_OUT_3*IFTIME>s CUMU_TT[COMP_3] THEN C_REL_3*1 0000 ELSE 0 
«# CR_INPUT_3 = if TIME > CUMU_TTICOMP_~3] then 0 else CR.RATE 

□ C_REL 4(t) = C_REL_4(t - dt) + (CR OUT_3 - CR.OUT.4 - CRJNPUT.4) * dt 
INIT C_REL_4 = 0 

INFLOWS: 

CR_OUT_3 = IF TIME >= CUMU_TT[COMP_3] THEN C_REL_3*10000 ELSE 0 
OUTFLOWS: 

CR_OUT_4 * IF TIME CUMU_TT[COMP_4] THEN C_REL_4*1 0000 ELSE 0 
^ CR JNPUT_4 = if TIME>CUMU_TT[COMP_4] then 0 else CR.RATE 

□ C_REL_5(t) = C_REL_5(t - dt) + (CR_OUT_4 - CR_OUT_5 - CRJNPUT.5) * dt 
INIT C_REL_5 = 0 

INFLOWS: 

CR_OUT_4 * IF TIME >* CUMU_TT[COMP_4] THEN C_REL_4'10000 ELSE 0 
OUTFLOWS: 

^ CR_OUT_5 * IF TIME > s CUMU_TT[COMP_5] THEN C_REL_5*1 0000 ELSE 0 
^ CR_INPUT_5 = if TIME>CUMUJT[COMP_5] then 0 else CR_RATE 

□ C_REL 6(t) = C_REL 6(t - dt) + (CR OUT_5) * dt 
INIT C_REL_6 = 0 

INFLOWS: 

^ CR OUT_5 * IF TIME >* CUMU_TT[COMP_5] THEN C_REL_5*10000 ELSE 0 



1 13 



I_ MASS_1(t) = MASS_1(t-dt) + (CR INPUT 1 - MASS OUT 1 - DISS PRECIP 1>- 
INIT MASS_1 = DOSE " " " - ' 

INFLOWS: 

CR_INPUT_1 = if TIME>CUMU_TT[COMP 1] then 0 else CR RATE 
OUTFLOWS: 

«$> MASS_OUT_1 = MASS_1*TRANSFERS[COMP 1] 
DISS_PRECIP_1 = ADJ_DISS_PRECIP[COMPj] 

□ MASS_2(t) = MASS_2(t-dt) + (MASS_OUT 1 + CR INPUT 2 • MASS OUT 2- 
DISS_PRECIP_2) • dt - - 
INIT MASS_2 = 0 

INFLOWS: 

MASS_OUT_1 = MASS_1*TRANSFERS[COMP_1] 
CRJNPUT_2 = if TIME>CUMU_TT[COMP_2] then 0 else CR RATE 
OUTFLOWS: 

MASS_OUT_2 = MASS_2*TRANSFERS[COMP_2] 
DISS_PRECIP_2 = ADJ_DISS_PRECIP[C0MP_2] 

□ MASS_3(t) = MASS_3(t-dt) + (CRJNPUT 3 ♦ MASS OUT 2 -MASS OUT 3- 
DISS_PRECIP_3) * dt ~ " " 

INIT MASS.3 = 0 
INFLOWS: 

CR_INPUT_3 = ifTIME>CUMU TT[COMP_3] then 0 else CR_RATE 
MASSJDUT.2 = MASS_2TRANSFERS[COMP 21 

OUTFLOWS: 

# MASS_OUT_3 = MASS_3*TRANSFERS[COMP_3] 
«J> DISS_PRECIP_3*ADJ_DlSS_PRECIP[COMP 3] 

□ MASS_4(t) = MASS_4(t - dt) + (CR_INP(JT_4 + MASS OUT 3 - MASS OUT 4 - 
DfSS_PRECIP 4)*dt " ~ 

init mass 4 = 0 
inflows" 

CR_INPUT_4 a if TIME>CUMU_TT[C0MP_4] then 0 else CR.RATE 
MASS_OUT_3 = MASS_3TRANSFERS[COMP 3] 

OUTFLOWS: 

# MASS_OUT_4 = MASS_4TRANSFERS{COMP 4] 
DISS_PRECIP_4 = ADJ_DISS_PRECIP[COMPj4] 

□ MASS_5(t) = MASS_5(t - dt) + (CRJNPUT 5 + MASS OUT 4 - MASS OUT 5 - 
DISS_PRECIP 5)*dt ~ ~ " 

INIT MASS_5 = 0 
INFLOWS: 

CR_INPUT_5 = if TIME>CUMU_TT[C0MP_5] then 0 else CR_RATE 

# MASS_OUT_4 = MASS 4*TRANSFERS[COMP 4] 

OUTFLOWS: 

MASS_OUT_5* if time>4 then MASS 5*TRANSFERS[COMP 5] else 0 
DISS_PRECIP_5 = ADJ_DISS_PREOP[COMP 5] 

□ MASS_6(t) = MASS 6(t - dt) + (MASS OUT 5) • dt 

init mass 6 = 0 
inflows" 
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MASS_0UT_5 = if time>4 then MASS_5*TRANSFERSIC0MP_5] else 0 
— ] S0L_1(t) = S0L_1(t - dt) + (DISS_PRECIP_1 - S0L_0UT_1 - FLUXJ) 4 dt 
I NIT S0L_1 = 0 
INFLOWS: 

DISS_PRECIPJ = ADJ_DISS_PRECIP[COMP_1] 
OUTFLOWS: 
^ SOL_OUT_1 = SOL_1'TRANSFERS[COMP_1] 
FLUX_1 = 

CONCENTRATIONS(COMP_1]*ADJ_PERM[COMPjrSURFACE_AREA[COMP_1] 

□ SOL_2(t) = SOL_2(t - dt) + (SOLJDUTJ + DISS_PRECIP_2 - SOL_OUT_2 - FLUX_2) • dt 
INITSOL_2 = 0 

INFLOWS: 

SOL_OUT_1 = SOL_1*TRANSFERS[COMP_1] 

DISS_PRECIP_2 = ADJ_DISS_PRECIP[COMP_2] 
OUTFLOWS: 
1$ SOL_OUT_2 = SOL_2*TRANSFERS(COMP_2] 
^ FLUX 2 = 

CONCENTRATIONS[COMP_2rADJ_PERMtCOMP_2rSURFACE_AREA[COMP_21 

□ SOL_3(t) = SOL_3(t - dt) + (DISS_PRECIP_3 + SOL_OUT_2 - SOL_OUT_3 - FLUX_3) * dt 
INIT SOL_3 = 0 

INFLOWS: 

DISS_PRECIP_3 = ADJ_DISS_PRECIP[COMP_3] 
^ SOL_OUT_2 ■ SOL_2TRANSFERS[COMP_2] 
OUTFLOWS: 
<f SOL_OUT_3 = SOL_3TRANSFERS[COMP_3] 

FLUX_3 = 

CONCENTRATIONS[COMP_3]*ADJ_PERM[COMP_31*SURFACE_AREA(COMP_3] 

□ SOL_4(t) = SOL_4(t • dt) + (DISS_PRECIP_4 + SOL_OUT_3 - SOL_OUT_4 - FLUX_4) * dt 
INIT SOL_4 s o 

INFLOWS: 

<g DISS_PRECIP_4 = AOJ_DISS_PRECIP[COMP_4] 
■Jf SOL_OLIT_3 = SOL_3*TRANSFERS{COMP_3] 
OUTFLOWS: " 

SOL_OUT_4 = SOL_4TRANSFERS{COMP_4] 
•» FLUX_4=" 

CONCENTRATIONS[COMP_4rADJ_PERMtCOMP_4rSURFACE_AREA[COMP_4] 

C SOL_5(t) = SOL_5(t - dt) + (DISS_PRECIP_5 + SOL_OUT4 - SOL_OUT_5 - FLUX_5) * dt 
INIT SOL_5 = 0 

INFLOWS: 

US 



DISS_PRECIP_5 * ADJ_DISS_PRECIP[C0MP_5] 
S0L_0UT_4 =~SOL_4*TRANSFERS[COMP_4] 

OUTFLOWS: 

SOL_OUT_5 = if time>4 then SOL_5*TRANSFERS(COMP_5] else 0 

# FLUX_5 = if time<32 then 
CONCENTRATIONS(COMP_5]*ADJ_PERM[COMP_5]*SURFACE_AREA[COMP_5]*(32-ti 
me)/48*(VOL_5/1 7.2) else 0 

C SOL_6(t) = SOL_6(t - dt) ♦ (SOL_OUT_5) * dt 
I NIT SOL_6 = 0 
INFLOWS: 

SOL_OUT_5 = if time>4 then SOL_5*TRANSFERS[COMP_5] else 0 
Q VOL_1 (t) = VOL_1 (t - dt) + (- REABSJ - VOL_OUT_1 ) * dt 
INIT VOL_1 = INlT_VOLUME 
OUTFLOWS: 
<^> REABS.1 =VOL_rVOL_PARM[COMP_1] 
VOL_OUT_1 =VOL^1*TRANSFERS[COMP_1] 
VOL_2(t) = VOL_2(t - dt) ♦ (VOL_OUT_1 - VOL_OUT_2 - REABS_2) * dt 
INITVOL_2 = 0 
INFLOWS: 

<g VOL_OUT_1 =VOL_rTRANSFERS[COMP_1] 
OUTFLOWS: 
<$> VOL_OUT_2 = VOL_2*TRANSFERS[COMP_2] 

# REABS.2 = VOL_2*VOL_PARM[COMP_2] 

Q VOL_3(t) = VOL_3(t - dt) + (VOL_OUT_2 - VOL_OUT_3 - REABS_3) * dt 
. INIT VOL_3 = 0 
INFLOWS: 

<g VOL_OUT_2 = VOL_2*TRANSFERS[COMP_2] 
OUTFLOWS: 
<g VOL_OUT_3 = VOL_3*TRANSFERS{COMP_3] 
REABS_3 = VOL_3*VOL_PARMtCOMP_3] 

□ VOL_4(t) = VOL_4(t - dt) + (VOl_OUT_3 - VOL_OUT_4 - REABS.4) • dt 
INITVOL_4 = 0" 

INFLOWS: 

•gp VOL_OUT_3 = VOL_3TRANSFERS[COMP_3] 
OUTFLOWS: 

VOL_OUT_4 = VOL_4*TRANSFERS[COMP_4] 
■£> REABS_4 « VOL_4*VOL_PARM[COMP_4] 

□ VOL_5(t) = VOL_5{t - dt) ♦ (VOL_OUT_4 • VOL_OUT_5 - REABS_5) * dt 
INIT VOL_5 = 0~ 

INFLOWS: 

•jf VOL_OUT_4 = VOL_4*TRANSFERS|COMP_4] 
OUTFLOWS: 

VOL_OUT_5 = VOL_5*TRANSFERS[COMP_5] 
REABS_5 = VOL_5*VOL_PARMtCOMP_5] 
[3 VOL_6(t) = VOL_6(t - dt) + (VOL_OUT_5) * dt 
INIT VOL_6 = 0" 



INFLOWS: 

# VOL_OUT_5 - VOL_5*TRANSFERS[COMP 5] 

— v OL_ABS_1(t) = VOL_ABS_1(t-dt) + (REABS 1) "dt 
I NIT VOL_ABS_1 =0 

INFLOWS: 

«*> REABSJ =V0L_rV0L_PARM[C0MP_1] 

VOL_ABS_2(t) = VOL_ABS_2(t-dt) + (REABS 2) * dt 
INITVOL ABS 2 = 0 

INFLOWS: 

# REABS.2 = VOL_2*VOL_PARM[COMP_2] 

Q VOL_ABS_3(t) = VOL_ABS_3(t-dt) + (REABS 3)*dt 
INIT VOL_ABS_3 = 0 

INFLOWS: 

REABS_3 = VOL_3*VOL PARM[COMP 3] 

□ VOL_ABS_4(t)sVOL_ABS_4(t-dt) + (REABS*4)'dt 
INITVOL ABS_4 = 0 
_ INFLOWS: 

% REABS_4 = VOL_4*VOL_PARM[COMP_4] 

£ C v °L-ABS_5(t) = VOL_ABS_5(t-dt) + (REABS 5)*dt 
m INITVOL_ABS_5 = 0 " 

INFLOWS: 

jjj # REABS_5 = VOL_5*VOL PARMfCOMP 5] 

yg MULTI DOSE CALCULATION 
yfij OUTPUT CALCULATIONS 

= O CR_Release(t) = CR_Release(t - dt) + (CR cumrate) * dt 
M INIT CR_Release = 0 
H INFLOWS: 

jU * CR_cumrate * CRJNPUT_1 + CRJNPUT_2+CR_INPUT_3 + CRJNPUT_4 + CR_INPUT_5 

g O CU^_D^ 

INFLOWS: 

CUMM_DISS RATE = 

DISS_PRECIP_1 + DISS_PRECIP_2 + DISS.PRECIP_3 + DISS_PRECIP_4 + DISS.PRECIP 

O ABSORBED.TOTAL = ABSORPTION>ABSORPTION_3 + ABSORPTION_4MBSORPTION_5 

O F Dp%=ABSORBED_TOTAUDOSE*100 
O FLUX_TOTAL a 'FLUX_2+FLUX 3+FLUX 4+FLUX 5 
g) PARMS ~ 

Q DOSE = 1000 
O INIT_VOLUME = 100 
O PARACELLULAR = 1 
O PH[C0MP_1] = 1.5 
O pH[COMP~2] = 5 
O pH[COMP_3] * 6.5 
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2 pH[C0MP_4] = 7 
J pH[C0MP_5] = 6.5 

^ SURF ACE_AREA[C0MP_1] = if PARACELLULAR =0 then 50*SA_ADJ[COMP 1]else 
50*SA_ADJ[COMP_1] 

0 SURFACE_AREA[COMP_2] = if PARACELLULARsQ then 1 50*SA_ADJ[COMP_2] else 
~ 7.5*SA_ADJ[COMP_2] 

0 SURFACE_AREA[COMP_3] = if PARACELLULAR=Q then 1000*SA_ADJ{COMP_3] else 
50*SA_ADJ(COMP_3] 

0 SURFACE_AREA[COMP_4] = if PARACELLULAR=0 then 1000*SA_ADJICOMP 4] eise 
50*SA_ADJ{COMP_4] 

0 SURFACE_AREA[COMP_5] = if PARACELLULAR=0 then 850*SA_ADJ[COMP_5] eise 

W 42.5 , SA_ADJ[COMP_5] 

0 TlME_IN_HOURS * TIME 

Q VOLPARMfCOMPJ] = 0*VOL_ADJ[COMPJ] 

0 VOL~PARM[COMP_2] = 0*VOL_ADJ[COMP_2] 
O 0 VOL_PARM(COMP_3] = 1.75*VOL_ADJ[COMP_3] 
« 0 VOL_PARM(COMP_4]« 1.75*VOL_ADJ[COMP_4] 
£ 0 VOlTPARM[COMP_5] = 0.10*VOL_ADJ[COMP_5] 
fe PERMEABILITY CALCULATION 
~ 0 ACT PE[COMPS] = [0 , 
*\ ' 0. " 

fi o. 

0, 

L O ADJ_PERM[COMP_1] = 

m (2/(1+EFFLUX AdJ[COMP 1]))*REGIONAL[COMP 1]*FLUX ADJ[COMP 1]'3600+(CARRIER_ 
LI DJ(COMP_1 ]*ACT_PElCOMP_1 ]*3600/(1 +(CONCENTRATIONS[COMP_lT/(Km[COMP_1 ] ))))"0 

S O ADJ_PERM(COMP_2] * 

(2/{1+EFFLUX_ADJtCOMP_2]))*REGIONAL[COMP 2]*FLUX ADJfCOMP 2]*3600+(CARRIER_ 
DJ[COMP_2rACT_PE(COMP_2r3600/(1+(CONCENTRATIONStCOMP_2]/(Km[COMP_2])))) 

0 ADJ_PERM[COMP_3] * 

(2/(1+EFFLUX ADJ(COMP 3]))'REGIONALlCOMP 3rFLUX_ADJ[COMP_3]'3600+(CARRIER_ 
DJ[COMP_3rACT_PEtCOMP_3]*3600/(1+(CONCENTRATIONSlCOMP_3l/(Km[COMP_3])))) 

O ADJ PERM(COMP 4] * 

(2/(1 +EFFLUX ADJICOMP 4]))*REGIONAL(COMP_4]*FLUX_ADJ{COMP_4r3600+(CARRIER_ 
DJ[COMP_4]*ACT_PE[COMP>] , 3600/(1+(CONCENTP^TIONS[COMP_4]/(Krn(COMP_4])))) 

0 ADJ PERM[COMP 5] ~ 

(2/(1+EFFLUX ADJ{COMP_5]))*REGIONAL[COMP_5]*FLUX_ADJ(COMP_5]*3600+(CARRIER_ 
DJ[COMP_5]*ACT PEfCOMP 5]*3600/(1 +(CONCENTRATIONS[COMP_5}/(Km[COMP_5])))) 



C Km[COMPS] * [1 , 

1 , 
1 , 

1 . 

1] 

. . PASS_PE(COMPS] s [0 , 
1.10E-06, 
2.17E-06 , 
4.06E-06 , 
3.80E-06 ] 

O RC[COMP_1] « PASS_PE[COMP 1]*0 



C «C[C0MP_21 = IF PASS_PE[COMP 2]>0 THEN 1 ELSE 0 
O RC[COMP_3] = IF PASS.PE[COMPl3]>0 THEN 2 ELSE 0 
O RC[COMP_4] = IF PASS_PE[COMP 4]>0 THEN 4 ELSE 0 
C RC[COMP_5] = PASS_PE[COMP 5r0 
O RCSUM = RC[COMP_2]+RC[COMP 3]+RC[COMP 4] 

O REGIONAL[COMP_1] = PASS_PE[COMP_1] + RCSUM*0 
O REGIONAL[COMP_2J * if RCSUM=2 then 

!f " 9 -° 1 1 926 * 2 594378 *LOGN(1/PASS PE[COMP 21) -0 06551 5 

*LOGN(1/PASS_PE[COMP 2D*2))*(-1 ) else ~ 

ifRCSUM=4then 

ifRCSUM=6then ~ 

PASS_PE[COMP 2] " 
O REGlONAL[COMP_3] « if RCSUM*1 then 

(EXP( -3.238469 + 1.509131 *LOGN(1/PASS PEfCOMP 21) -0 022109 
•LOGN(1/PASS_PE[COMP 2)^2)^) else " tl ° OMP - 2]) -° 022109 
if RCSUM=4 then 

(EXP(-0.093739'LOGN(1/PASS_PE[COMP 4])+0.182344«LOGNn/PASS PFrroMP am n 
23631-LOGN(1/PASS PE{COMP 4 Wm) etee LOGN ( 1 /PASS_PE[COMP_4])A 2 -0.00 

ifRCSUM=5then ~ 
PASS_PE[COMP_3] 



r, REGI0NAL(C0MP_4] * if RCSUM-1 then 

<EXP( 14.455255 -1 .264630 *LOGN(1/PASS_PE[COMP 2]) + 0.06201 5 
•LOGN(1/PASS_PE[COMP_2]) A 2)) A (-1 ) else 
if RCSUM=2 then 

(EXP( 1 1 .480333 -0.791 109 *LOGN(1/PASS_PE(COMP_3l) + 0.066063 
*LOGN(1/PASS_PE[COMP_3]) A 2)) A (-1 ) else 
if RCSUM=3 then 

0.5*(EXP( 14.455255 -1 .264630 *LOGN(1/PASS_PE[COMP_2]) ♦ 0.082015 
*LOGN(1/PASS_PE[COMP_2]) A 2)) A (-1) 

+0.5*(EXP( 1 1 .480333 -0.791 109 *LOGN(1/PASS PE[COMP_3]) ♦ 0.066063 

•LOGNll/PASS.PElCOMP.SJ^jr^l ) else 

PASS PEfCOMP 4] 
Q REGIONAL[COMP_5] - PASS_PE[GOMP_5] +RCSUM*0 
O) SOLUBILIY CALCULATION 

Q ADJ_SOLUB[COMP_1] - if USER_pH[COMP_1]>spH[COMP_1} then USER_SOLUB[COMP_1] 

else 

ri ((USER_SOLUBlCOMP 2]-USER SOLUB[COMP_1])/(USER_pH{COMP 2]-USER_pH(COMP 1] 

))'(pH[COMP 1]-USER_pH[COMP_1])+USER_SOLUB(COMP_1] 
.n O ADJ SOLUB[COMP_2]»ifUSER_pH[COMP_2]*pH[COMP_2]thenUSER_SOLUB[COMP_2l 
m else if USER pH[COMP_2]<pH[COMP_2] then 

J ((USER_SOLUB[COMP 3J-USER SOLUB[COMP_2])/(USER_pH[COMP 3RISER pH[COMP 2] 
§ ))*(pH[COMP 2]-USER_pH[COMP 2])+USER_SOLUB{COMP_2} else 

Qi ((USER_SOLUB[COMP 2]-USER SOLUB[COMPJMUSER_pH[COMP 2]-USER pH[C0MP 1} 
UJ )r(pH[COMP_2]-USERj)HlCOMP_1])+USER_SOLUB{COMP_1J 

7 O ADJ_SOLUB[COMP_3] = if USER_j>H[COMP_3]*pH[COMP_3] then USER_SOLUB(COMP 3] 
u else if USERj>HtCOMP_3]<pH[COMP_3J then 

U ((USER SOLUBfCOMP 4].USER_SOLUB[COMP_3])/(USER_pH[COMP_4].USER_pH[COMP 3] 
fU )r(pH[COMP_3]-USERj)H{COMP_3])+USER_SOluB[COMP_3] else 

U ((USER SOLUB[COMP 3]-USER SOLUB[CO"MP 2])/(USER_pH[COMP_3]-USER pH[COMP 2] 
O ))*(pH[COMP_3]-USER_pHtCOMP_2])+USER_SOLUB{COMP_2] 

M> Q ADJ SOLUBICOMP 4]aifUSER_pH[COMP 4]apH{COMP_4]thenUSER_SOLUB[COMP_41 
elselfUSER pHfCOMP 4]<pH[COMP 4] then 

((USER_SOLUB[COMP 5J-USER SOLUB[COMP_4])/(USER _pHlCOMP_5]-USER_pH[COMP_4] 
))*(pH(COMP_4]-USERj>H[COMP 4])+USER SOLUB(COMP_4] else 

((USER SOLUB[COMP 4J-USER SOLUBFCOMP_3])/(USER _pH[COMP_4]-USER_pH[COMP_3] 
)r(pH[COMP 4].USER3HtCOMP_3])+USER_SOLUB[COMP_3] 
Q ADJ_SOLUB[COMP_53 * if USER_pHlCOMP_3]*pHlCOMPJ3] then USER_SOLUB{COMP_3] 
elself USER_pH[COMP_3]<pH{COMP_3] then 

{(USER SOLUB[COMP 4J-USER SOLUB[COMP_3]y(USER_pH[COMP_4]-USER_pH[COMP_3] 
))*(pHtCOMP 3]-USER_pH[COMP 3])+USER_SOLUB[COMP_3] else 

((USER SOLUBICOMP 3]-USER SOLUB(COMP_2]y(USER_pH[COMP_3]-USER_pHtCOMP_2] 
))*(pH[C~OMP 3J-USER pH(COMP_2])+USER_SOLUB[COMP_2] 
O USER_pH[COMPS] = [1 .5 , 
5. 

6.5, 

7. 

7.5] 



\L0 



~, USER_SOLUB(COMPS] = [31 , 

3.65," 
3.65 . 
3.65 . 
3.65] 
>□) TRANSIT TIME 

Q ADJ_TRANSIT_TIME[C0MP_1] = .5*TRANSIT_ADJ[C0MP_1]*USER_TT_INPUT 
J ADJ_TRANSIT_TIME[C0MP_2] = .25*TRANSIT_ADJ[COMP_2]*USER_TT_INPUT 
"j ADJ_TRANSIT_TIME[C0MP_3] = 1.5*TRANSIT_ADJ[COMP_3rUSERlTTJNPUT 
(f ADJ_TRANSIT_TIMEIC0MP_4] * 1 .5*TRANSIT_ADJ[COMP_4]*USER_TTJNPUT 
Q ADJ J"RANSIT_TIME[C0MP_5] = 24'TRANSIT_ADJ[COMP_5]«USER_TT_INPUT 
Q CUMU_TT[C0MP_1] = ADJ_TRANSIT_TIME[C0MP_1] 

Q CUMU_TT[COMP_2] = ADJ_TRANSIT_TIMEICOMP_1]+ADJ_TRANSIT_TIME[COMP_2] 
O 0 CUMU_TT[COMP 3] = 

m ADJ TRANSIT TlME{COMP 1J+ADJ TRANSIT TIME[COMP 2J+ADJ TRANSIT_TIME[COMP 

£ 3] 

08 Q CUMU TT{COMP 4] = 

m ADJ T^NSIT_TIMEIC0MP_1]+ADJ_TRANSIT^TIME[C0MP_2]+ADJ_TRANSIT_TIME[C0MP 
tn 3]MDJ_TRANSITTIMEIC0MP_4] 
W Q CUMU_TT[COMP 5] = 

«J ADJ_TRANSIT TIMEfCOMP 1]+ADJJ"RANSITTIME[COMP_2]+ADJ_TRANSITTIME[COMP 
3]+ADJ_TRANSIT TIME[COMP 4]+ADJ_TRANSIT_TIME[COMP_5] 
O TRANSFERSICOMPJ] = LOGN(10)/ADJ_TRANSlf_TIME[COMP_1] 
tf. O TRANSFERSfCOMP.2] - LOGN(10)/ADJ_TRANSIT_TIME[COMP_2] 
Jj O TRANSFERS{COMP_3] = LOGN(10)/ADJ_TRANSITTIME(COMP_3] 
S Q TRANSFERS[COMP~4] = LOGN(10)/ADJ_TRANSITTIME[COMP_4] 
S O TRANSFERS[COMP_5] * LOGN(10)/ADJ_TRANSIT_TIME(COMP_5] 
*~ Q USERJTJNPUT = 7 
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